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ارزبایی کارایی روش پیش‌بینی موضعی ناشی از بازسازی فضای فاز و مدل شبکاً عصبی مصنوعی به‌منظور 
پیش‌بینی رواناب (مورد مطالعاتی: ایستگاه پل کهنه» کرمانشاه)* 


ی ۰ (۱) . ۳ 
محمد ذونعمت کرمانی خاطره امیرخانی 


چکیده در اين تحفیق دور روش شبکة عصبی مصنوعی و روش پیش‌بینی موضعی بر پایةٍ تلور ی آشوب به‌منظور پیش‌بینی روزانه, هفنگی و 
ماهانةٌ رواناب ایستگاه پ لکهنه بر رودخحانهة قره‌سو به‌کا رگرفته شده‌اند. نتایج حاصل از پیش‌بینی با استفاده از روش پیث‌بینی موضعی مبین 
نزدیکی نتایج با داده‌های مشاهداتی در مقیاس‌های روزانه و هفتگی و عدم تطابق مناسب در مقیاس ماهانه بوده است که بیانگر وجود 
آشوبناکی در مقیاس‌های روزانه و هفتگی است. هم‌چنین تنایج پیثربینی شده با استفاده از شبکا عصبی در مقیاس‌های روزانه و هفتکی 
برحلاف مقیاس ماهانه ضعیف‌تر از روش پیش‌بینی موضعی بوده است. 


واژه‌های کلیدی سری زمانی رواناب نظری؛ آشوب پیش‌بینی موضعی, شبکة عصبی مصنوعی. 
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* تاریخ دریافت مقاله ٩۳/۸/۲۱‏ و تاریخ پذیرش آن ۹۶/۸/۲۶ می‌باشد. 
(۱) نویسنده مسئول دانشیار بخش مهندسی آب. دانشگاه شهید باهنر کرمان. )ما20 


(۲) دانش آموخحتهة کارشناسی ارشد. مهندسی منابع آب دانشگاه شهید باهنر کرمان. 


۹۲ 


مقد مه 
یکی از مهم‌ترین مطالعاتی که در هیدرولوژی انجام 
می‌شود. بررسی و تعیین مقدار رواناب و برآورد حداکثر 
سیلاب و تهية آب‌نمود آن می‌باشد. هم‌چنین برآورد 
رواناب حاصل از رواناب‌های جوی از اهمیت خاصی 
در هیدرولوژی و مدیریت حوضه‌های آبخیز و حفاظت 
از آب و خاک برخوردار است. مسائل مربوط به 
رواناب. گوناگون است و طبیعت بسیار پیچیده‌ای دارد. 
به‌دلیل وقوع رواناب در حوض آبریز رودخانه‌ها 
فرسایش‌های عظیمی در سطح حوضه اتفاق می‌افتد و 
نهایتاً مشکلات عدیده‌ای را در پایین‌دست به‌وجود 
می‌آورد که ناشی از تهنشینی و انباشته شندن رشوب 
است. برای کنترل و هدایت رواناب و تخلیةُ جریان‌های 
سطحی در اراضی کشاورزی و انتقال آنهابه محل 
مناسب و خارج از منطقه مدل‌های متفاوتی توسط 
محققان. سازمان‌های مطالعاتی- تحقیقاتی در کشورهای 
مختلف جهان ارائه شده و مورد استفاده نیز قرار گرفته 
است. از آن جمله می‌توان به روش‌های تجربی و 
غیرمستقيم اشاره کرد که هر کدام دارای یک‌سری 
ضرایب منطقه‌ای خاص به‌خود بوده و تحت آن شرایط 
توسعه داده شده‌اند. به اين دلیل در سایر مناطق ممکن 
است غیرقابل استفاده باشد یا استفاده از آنها همراه با 
خطاهایی باشد. باوجود روابط غیرخطی, عدم قطعیت و 
عدم صراحت زیاد و ویژگی‌های متغیر زمانی و مکانی 
در سیستم‌های گردش آبی. هیچ‌یک از مدل‌های آماری و 
مفهومی پيشنهادشده به‌منظور الگوسازی دقیق نتوانسته 
به‌عنوان یک مدل برتر و توانا شناخته شوند. امروزه 
شبکه‌های غیر خطی به‌عنوان یکی از سیستم‌های 
هوشمند در پیش‌بینی یک چنین پدیده‌های پیچیده بسیار 
مورد استفاده قرار می‌گیرند [۱]. از آن‌جمله شبکه‌های 
عصبی مصنوعی به‌عنوان مدل‌های جعبه‌سیاه ابزاری 
توانمند برای مدل‌سازی سیستم‌های غیرخطی هستند و 
قادرند روابط پیچیدهُ حاکم بر ایین فرآیندها را مدل 


مان زدیا توت کرسای و همکاران(۲۹۱۲۳) اه 


نشريهُ مهندسی عمران فردوسی 


ارزیاب ی کارایی روش پیش‌بینی موضعی ناشی از بازسازی فضای فاز و ... 


پیگن‌یتی زواناتب روزاه با اشفاهه از ووش‌هان هه 
عصبی مصنوعی پیشرو و تابع شعاعی پرداختند [3] 
نتایج حاصل از تحقیق برتری روش غیرخطی شبکه 
عصبی را نسبت به روش‌های خطی محرز می‌سازد. یکی 
دیگر از شیوه‌های مدل‌سازی غیرخحطی, استفاده از 
مفاهیم تئوری آشوب می‌باشد. نظري؛ آشوب به‌دلیل 
داشتن خاصیت دینامیکی غیرخطی و قطعی بودن مورد 
توجه خاص در بررسی فرآیندهای هیدرولوژیکی بوده 
است [۶] هم‌چنین پیش‌بینی موضعی که مبتنی بر 
پارامترهای نظریة آشوب می‌باشد یک روش جدید در 
پیش‌بینی چنین پارامترهایی است. آشوب در لغت 
به‌معنای هرج و مرج و بی‌نظمی است. انگار اصلی و 
کلیدی تئوری آشوب این است که هر بی‌نظمی. نظمی 
نهفته است. به این معنا که نباید نظم را تنها در یک 
مقیاس جستجو کرد پدیده‌ای که در مقیاس محلی. 
کاملاً تصادفی و غیر قابل پیش‌بینی به‌نظر می‌رسد چه‌بسا 
در مقیاس بزرگتی کاملاً پایا و قابل پیش‌بینی باشد. سه 
مزیت عمده از نظرية آشوب برای تحلیل سیستم‌های 
هیدرولوژیکی وجود دارد: 
) برای سیستم شدیداً غیرخطی قابل اجراست. 
۲) به‌طور طبیعی کلب دینامیک‌های مهم سیستم را 
محاسبه می‌نماید. 
۴ اطلاعات و روابط سیستم را بدون نیاز به کشف 
قوانین يا روابط تحت دینامیک کشف می‌کند. اگرچه 
تحلیل آشوب. یک رابطه برای استخراج مدل متداول 
تهیه نمی کند. اما اطلاعات ارزشمندی در رابطه با همه 
سیستم‌های دینامیکی برای کاربرد کلی از قبیل پیش‌بینی 
داده‌ها استخراج می‌کند. نظريةٌ آشوب به‌نظر می‌رسد 
به‌طور طبیعی برای سیستم‌های هیدرولوژیکی قابل اجرا 
باشد [5] 

از آنجایی که پیش‌بینی و تعیین میزان کمی 
فرآیندهای تولید رواناب و انتقال آن به نقطه خروجی 
حوض آبخیز از اهمیت خاصی برخوردار است بررسی 


وضعیت دینامیکی سری زمانی رواناب می‌تواند در 


سال بیست و هشتم شمارهٌ دو ۱۳۹۵ 


محمد ذونعمت کرمانی- حاطره امیرحانی 


انتخاب و گزینش مدل‌های مناسب برای شبیه‌سازی و یا 
پر یی ان راکفا پاش هو همین زاستته سا پروسیتین 
تصادفی و یا آشوبناک بودن سری زمانی مورد نظر 
می‌توان به اطلاعات ارزشمندی در رابطه با نحوه رفتار 
این پدیده دست یافت و با استفاده از آن به بررسی 
روش‌های پیش‌بینی مختلف و یافتن روشی مناسب و 
دقیق به‌منظور پیش‌بینی رواناب دست یافت. ذیلا 
پژوهش‌ها و تحقیقات مرتبط در این زمینه معرفی و 
بررسی می‌گردند. 

جایاواردنا رواناب را براساس روش تقریب محلی 
مورد پیش‌بینی قرار داد [6]. نتایج حاصل این تحقیق 
نشان دادند که روش تقریب محلی می‌تواند روشی 
مناسب برای پیش‌بینی رواناب باشد در حالی‌که باید در 
مناطق دیگر که دارای توپوگرافی متفاوتی هستند این 
روش مورد آزمایش قرار گیرد. سیواکومار و همکاران 
آشسوبناکی فرآیند بارش -رواناب ماهانه در حوضة 
رودخانء گوتا در سوئد را بااستفاده از روش بعد 
همبستگی روی سری‌های بارش رواناب و ضرایب 
رواناب بررسی نمودند. با وجود رفتار آشوبی در سه 
نوع سری زمانی مورد مطالعه, نتیجه گرفتند که فرآیند 
بارش-رواناب حوض؛ مذکور رفتاری آشوبناک نشان 
می‌دهد [7]. در تحقیق دیگر سیواکومار و همکاران [8] 
به پیش‌بینی رواناب ماهانه در حوضهء رودخانهة 
اوراکیوری در برزیل با استفاده از بازسازی فضای فاز 
پرداختند. آنها از روش پیش‌بینی خحطی بااستفاده از 
چندجمله‌ای‌های محلی که بر اساس پارامترهای تئوری 
آشوب توسعه يافته. استفاده نمودند. نتایج نشان دادند که 
رواناب دارای رفتاری آشوبناک است و روش پیش‌بینی 
حطی باتوجه به ضریب همبستگی ۰/۸۹ و ضریب 
بهره‌وری ۰/۹۶ روشی مناسب برای پیش‌بینی رواناب 
می‌باشد [8]. ایسلام و سیواکومار [9]به بررسی 
خصوصیات و پیش‌بینی رواناب با یک دید دینامیکی 
غیر خطی پرداختند. آنها بااستفاده از داده‌های رواناب 


روزانة حوضهء لندربورگ در دانمارک ابتدا به بررسی 


سال بیست و هشتم. شماره دو ۱۳۹۵ 


٩۳ 


رفتار دینامیکی رواناب با روش‌های بعد همبستگی» 
طیف توانی تبدیل فوریه و بازسازی فضای فاز پرداختند 
و سپس با روش پیش‌بینی غیرخطی پیش‌بینی رواناب را 
انجام دادند. نتایج حاکی از وجود رفتار آشوبناک با بعد 
کم افرراین ری پزمانن ی میتی تسعا دقیقی با روش 
قر قظی برهد [ و راکو شنک را [10] سای 
خریال فاهات شا باماه و با امطاف رشان ری 
موضعی نظرية آشوب انجام دادند. آنها نشان دادند که 
برای پیش‌بینی‌های کوتاه‌مدت. این روش نتیجه بهتری 
نسبت به روش‌های دیگر داراست [10]. هو و هوانگ 
به بررسی آشوبناکی سری رواناب در ایستگاه فودیدین 
در چین پرداختند. بررسی‌های آنها حاکی از وجود 
نشانه‌های قانع کننده‌ای از آشوب در این سری زمانی 
راتس 1 

پژوهشگران کشورمان نیز تحقیقاتی را با استفاده از 
نظریة آشوب بر روی پدیده‌های مختلف هیدرولوژیکی 
انجام داده‌اند که از آن جمله می‌توان به مطالعات فرزین 
و همکاران درخصوص تحلیل آشوب‌پذیری رواناب 
مامانه [۱۲) قربانی و همکاران در رابطه با پیش‌بینی 
ذرات معلق رسوب در رودخانه [4] حسن‌زاده و 
همکاران در زمينة پایش خشکسالی [۱۳]» مرادی‌زاده و 
هیکازان در ارتاط با تس مرسان رویعاته ۶1 
انیس‌حسینی و ذاکر مشفق به‌منظور تحلیل فرآیند 
بارش -رواناب [۱۵] قربانی و همکاران درخصوص 
تحلیل رواناب [16] و جانی و همکاران در رابطه با 
تحیل خصوصیات دینامیکی بارش ماهانه [17] اشاره 
نمود. در این میان» انیس حسینی و ذاکر مشفق کاربرد 
نظریة آشوب در تحلیل فرآیند بارش-راوناب را در 
حوضء رودخانء کشکان بااستفاده از روش بعد 
همبستگی و توان لیاپانف مورد بررسی قرار دادند [۱۵] 
آنها از داده‌های روزانهٌ دبی رودخانهٌ کشکان بدین منظور 
استفاده نمودند. نتایج نشان دادند که شواهد قانع کننده‌ای 
از وجود آشوب در فرآیند رواناب- رواناب حوضٌ 


مذکور مشاهده می‌شود بنابراین شرایط برای پیش‌بینی 


نشريهُ مهندسی عمران فردوسی 


۹۶ 


رفتار این فرآیند به‌کمک روش‌های مبتنی بر نظرية 
آشوب مناسب تلقی می‌شود. هم‌چنین قربانی و همکاران 
کاربرد نظرية آشوب و الگوریتم ژنتیک در پیش‌بینی 
سری زمانی رواناب را مورد بررسی قرار دادند [116 
آنها نشان دادند که رواناب رفتاری آشوبناک دارد و 
پیش‌بینی این سری زمانی بااستفاده از الگوریتم ژنتیک از 
دقت بالایی برخوردار می‌باشد. 

با توجه به اهمیت پیش‌بینی سری زمانی رواناب و 
محدود بودن منابع تحقیقاتی مرتبط بااستفاده از نظرية 
آشوب در بررسی این پدیده در داحل کشور در مقالة 
حاضر به بررسی توانایی نظرب؛ آشوب در مدل‌سازی 
پدیدة رواناب پرداخته شد که در ابتدا آشوبناکی سری 
زمانی رواناب در مقیاس‌های زمانی روزانه» هفتگی و 
ماهانه در ایستگاه پل کهنه واقع در حوض آبریز قره‌سو 
بررسی گردید. سپس پیش‌بینی رواناب با بهره جستن از 
روش غیرخطی پیش‌بینی موضعی صورت پذیرفت. 
به‌منظور ارزیابی کارایی روش پیش‌بینی موضعی. 
پیش‌بینی رواناب بااستفاده از شبکه عصبی مصنوعی 
پرسپترون چندلایه نیز انجام شد و نتایج مورد بررسی 


قرار گرفت. 


ارزیاب ی کارایی روش پیش‌بینی موضعی ناشی از بازسازی فضای فاز و ... 


منطقةٌ مورد مطالعه 

حوض؛ آبریز قره‌سو به مختصات جغرافیایی ۳۳۰۰۳۲۵۳ 
تتا ۵۵۰۸ تفن الیو ۲۳۲۲۲ 39۱ 
۲۳ طول شرقی با مساحت ۵18۲/۷ کیلومتر 
مربع از حوضه‌های داخلی استان کرمانشاه می‌باشد. این 
حوضه از حوضه‌های سرشاخه کرخه است و از جنوب 
استان کردستان آغاز می‌شود و تا جنوب‌شرقی استان 
کرمانشاه ادامه می‌یابد [۱۸]. 

به‌منظور بررسی سری زمانی رواناب حوضه 
قره‌سی اطلاعات مورد نباز از سازمان آب منطقه‌ای 
غرب در استان کرمانشاه مربوط به دوره‌های آماری 
ایستگاه پل کهنه (داده‌های روزانه از تاریخ ۱۳۹۷/۰۷/۰۱ 
لغایت ۱۳۷۸/۰۳/۳۱) گرفته شد. مشخصات ایستگاه 
مورد نظر در جدول (۱) نشان داده شده است. هم چنین 
موقعیت جغرافیایی حوضه قره‌سو و ایستگاه پل کهنه در 
شکل (۱) قابل مشاهده می‌باشد. افزون بر این 
مشخصات آماری سری زمانی رواناب در ایستگاه 
پل کهنه در جدول (۲) و نمودار تغییرات رواناب در 
مقیاس ماهانه در شکل (۲) ارائه شده است. 


جدول ۱ مشخصات ایستگاه هیدرومتری پل کهنه (شرکت آب منطقه‌ای غرب. ۱۳۸۳) 


ایستگاه رودخانه 
طول جغرافیایی 


۶۷-۷ 


پل کهنه 


قره‌سو 


سال تأسیس 
عرض جغرافیایی ارتفاع از سطح دریا(متر) 
۳۶-۱ ۱۳۸۳ ۱۳۳۶ 


جدول ۲ مشخصات آماری سری زمانی رواناب در ایستگاه پل کهنه (کلیه واحدها بر مترمکعب بر ثانیه می‌باشد) 
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ایستگاه پل کهنه 
مقیاس زمانی 
میانگین | انحراف معیار | حداکثر مقدار | حداقل مقدار | تعداد داده 
و ۱ ۱۹/۲۰۷ ۳۳۸/۹۹ ۱ 1۳۹۲ 
هفتگی ۱۱۷ ۱۱۸ ۵ ۱۵ ۰ ۳۶ 
ماهانه ۹۳/۹۰ ۱۳:۸ ۳/۵۹( ۰ :۱ 
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۹6۵ 


شکل ۱ موقعیت حوض آبریز قره‌سو و ایستگاه هیدرومتری پل کهنه 


2500 - 

۱ 3 2000 ظً 

۱ 0 ؛ 

1 1000 - ۱ ۱ 

500 1۱ ۱ ۱ ۱ 

۱ ۱ ۱ ۱ ال مه در ال ۷ ۳ 

۱ / ار ال لا اه ۸۸ أله + 0 كِ 
150 100 50 0 


زمان (ماه) 


شکل ۲ نمودار سری زمانی رواناب در مقیاس ماهانه در ایستگاه پل کهنه در بازهُ زمانی سال‌های ۱۳۱۷ الی ۱۳۷۸ 


مواد و روش‌ها 
تحلیل غیرخطی بر پاية مفاهیم تئوری آشوب 
این بخش به بررسی و تحلیل غیرخطی براساس مفاهیم 
تلوری آشوب شامل بازسازی فضای فاز» تعیین بعد 
همبستگی و بزرگ‌ترین نمای لیاپانف می‌پردازد. باتوجه 
به نتایج حاصل از تحلیل موارد فوق می‌توان وضعیت 
آشوبناکی و یا تصادفی بودن سری زمانی را تعیین نمود. 


بازسازی فضای فاز. در بررسی دینامیک سیستم‌های 
آشوبناک» بازسازی فضای فاز فرآیندی مهم به‌شمار می- 
رود. روش‌های مختلفی ترا این مقوله وجود دارد که 
رایج‌ترین آنها برای بازسازی فضای فاز سری زمانی؛ 


روش ارائه‌شده توسط تاکنز می‌باشد [5]. در این روش 


سال بیست و هشتم. شماره دو ۱۳۹۵ 


ملفه‌های سیستم در یک فضای چندبعدی باتوجه به 
زمان تأخیر مناسب سیستم و بعد محاط سری زمانی 
نگاشت می‌شوند. به‌منظور تعیین زمان تأعیر می‌توان از 
روش‌های هم‌چون نمودار همبستگی‌نگار تابع اطلاعات 
میانگین متقایل و یا بعد همبستگی استفاده نمود که در 
این بین روش اطلاعات میانگین متقابل به‌علت توانایی 
در خذ خصوصیات غیرخطی سری زمانی بیشتر سورد 
توچه قرار می‌گیرد: افژون بر این جهت نمیین بمد بحاط 
می‌توان از روش تیک تر بر همسایگی کاذب استفاده 


نمود. 


الگوریتم بعد همبستگی گراسبرگ ر-پروکاشیا. آنالیز 
تابع همبستگی که به‌عنوان روش بعد همبستگی نیز 
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اه رد یکی از تکمک‌هاین اسب که رز 
گسترده برای بررسی وجود یا فقدان آشسوب در 
سری‌های زمانی استفاده می‌شود. از میان تعداد زیادی 
الگوریتم‌های موجود برای محاسبهٌ تابع همبستگی, 
الگوریتم گراسبرگر-پروکاشیا معمولاً بیش از سایر 
روش‌ها استفاده می‌شود [5]. الگوریتم گراسبرگر- 
پروکاشیا با استفاده از بازسازی فضای فاز سری زمانی 
مقدار بعد همبستگی را تعیین می‌نماید. فضای فاز یک 
فضای انتزاعی است که میزان درجه آزادی حرکت 
سیستم می‌باشد. برای سری زمانی اسکالر ,26 که در آن 
آالر...,1-1,2,3 فضای فاز را می‌توان با استفاده از روش 
زمان تأخیر به‌دست آورد [7]. در روش زمان تأخیر 1 
و زمان‌های بعدی آن با هم ترکیب می‌شوند و مختصات 
یک سری زمانی برداری را تعیین می‌کنند [19]: 

که در آن ۳0-1(/۸)- ...1,2,3 ۲ بعدی 
از بردار ۷ می‌باشد که بعد محاط نام دارد و ۲ به‌عنوان 
پارامتر زمان تأخیر بیان می‌شود. برای فضای فاز 17 
بعدی تابع همبستگی (065) به‌صورت زير می‌باشد: 


(اکزعکن 


)0 ۳-۳ 2 وسب زا -0* 


که در آن ۲1 تابع هویساید پله‌ای با 1-(1100 برای 
0 و ۲۲)۱(<0 رای 0 گ ۱۷1 مسئی‌باشد 
و ارلا- 2۲-۷ 0 ۰ آ[ تعداد تقاط در فضای مزبور 1 


شعاع کره به مرکز :۷ يا ز۷ می‌باشد. فاصلهٌ 


بین بردار تأخیر 0 بعدی به‌دست آمده از معادلة (۱) 


می‌باشد. چنان‌چه سری زمانی به‌وسیله یک جاذب 


توصیف شود پس انتگرال همبستگی (0 با شعاع ۲ 

به‌صورت زير وابسته می‌شود: 

)۳( ع< ور(00) 
سس۱۳ 


در این رابطه ۰2 یک مقدار ثابت است و ۲ توان 


همبستگی یا شیب منحنی (6) 108 در مقابل (1086 


می‌باشد: که تشأن‌دهنده میزان اتخزافت یک سرخ زمانی از 
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ارزیاب ی کارایی روش پیش‌بینی موضعی ناشی از بازسازی فضای فاز و 


حالت تصادفی می‌باشد و به‌صورت زیر به‌دست 0 


۳ 10800 9 
10801 مب هت ِِ 


به‌طور کلی شیب به‌وسیلهة برازش یک خط مستقیم 
با روش حداقل مربعات داخل یک ناحية معین به نام 
ناحیةٌ مقیاس گذاری برآورد می‌گردد. اگر توان همبستگی 


کلی سیستم دارای دینامیک آشوبناک می‌باشد. مقدار 
اشباع از توان هنکن به‌عنوان بعد تسیچ () از 


تفای قعر یف مر شورو [ ۲19,20 


بزرگ‌ترین نمای لیاپانف. نماهای لیاپانف. مقدار 
نماهای متوسط از واگرایی یا همگرایی مدارها در فضای 
فاز می‌باشند. برخی از محققان بیان کرده‌اند که وجود 
یک نمای لیاپانوف مثبت حاکی از حضور آشوب در 
سیستم می‌باشد [21] درحالی که برخی دیگر بیان 
کرده‌اند که این شرط لازم است اما کافی نیست [22]. در 
اغلب موارد. یک سیستم با حداقل یک نمای لیاپانوف 
مثبت. آشوبناک درنظر گرفته می‌شود [23]. روزنستین 
(۱۹۹۳)» شیوة جدیدی برای محاسبهٌ بزرگ‌ترین نمای 
لیاپانوف بااستفاده از سری‌های زمانی ارائه نمود. این 
روش. بعد از بازسازی فضای فاز مناسب با استفاده از 
مقادیر 7و 0 یک نقطُ ,, را انتخاب و تمام نقاط 
همسایه «۷ با فاصله‌ای نزدیک‌تر از یک فاصلهٌ 7 از آن 
را پیدا می‌کند و متوسط فاصلهٌ آنها از آن نقطه را 
محاسبه می‌کند [24]. این روند برای محاسبةٌ یک مقدار 
متوسط 5 به نام فاکتور کشش. برای لا نقطه در طول 


1 


ره 


نقطةهً 5 می‌باشد. در سیستم‌های دارای رفتار بی‌نظم, 


۲ -< ]۰ بعنی نمودار ٩‏ بر حسب لا یا] نشان 


1 1 
دس 


تا تعداد همسایه‌های یافت‌شده در اطراف 
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خواهد داد که این منحنی در ابتدا به‌صورت خطی آغاز 
می‌شود و دارای یک ناحيةٌ مسطح در امتداد آن می‌باشد. 
شیب این نمودار نشان‌دهنده تقریب بزرگ‌ترین نمای 
لیاپانوف می‌باشد هم‌چنین باتوجه به بزرگ‌ترین نمای 
لیاپانوف به‌دست آمده مقدار افق پیش‌بینی (۸) 
به‌صورت زير محاسبه می‌شود: 

۸ - 2 62 


بر | نج 


است. در این تحقیق از روش روزنستین به‌منظور 
تتتاشته زر گتییی تما ابا انشفاهه شوه انیم 


فرایند پیش‌بینی 
در این تحقیق از دو روش غیرخحطی برای پیش‌بینی 
روانات استفاده شده است. روش بخست (نتسن نی 
موضعی) بر پا ایجاد فضای فاز و توجه به زمان تأخیر 
بنا شده استت! قی معالی .که در روش دوم که می‌توان 
یکی از رایج‌ترین روش‌های پیش‌بینی غیرخطی برشمرد. 
از شبکهة عصبی مصنوعی پرسپترون جندلایه استفاده 


شده است. 


پیش‌بینی موضعی. فرض می‌شود که داده‌های سری 
زمانی از یک سیستم دینامیک آشوبناک با بعد همبستگی 
۵ در فضای فاز تشکیل شده باشد. فضای فاز جدیدی را 
می‌توان با اضافه نمودن مختصات مستقل اضافی تشکیل 
داد. یکی از این مختصات‌های مستقل خود سری زمانی 
است و بقیةٌ مختصات‌ها با تأخیر دادن (0-1) برابر زمان 
تأخیر (یعنی (۳-1)) سری زمانی اصلی به‌دست 
می‌آیند به‌طوری‌که همبستگی بین مختصات‌ها صفر 
شود. در این صورت سری زمانی به‌صورت زیر درنظر 
گرفته می‌شود (کوجاک و همکاران» ۲۰۰۷): 

0( ار ان 


بازسازی به‌صورت زیر انجام می‌گیرد: 
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۹۷ 


1 یی فد ۳9 < ,6 
۷ با و,1(7-) + 1(,2-) +1 - ز 

(۷) 

که در آن :26 یک بردار 70 بعدی می‌باشد. زمانی که 
1 < 7 است تبدیل ذکرشده یک محاط‌سازی است 
که در آن فضای حالت کاذب ایجادشده. ساختار جاذب 
محاط‌شده در فضای حالت اصلی را حفظ می‌کند. در 
چنین فضای بعدی. پیش‌بینی به‌وسبله تخمین تغییر زک 
با زمان صورت می‌گيرد. با درنظر گرفتن ارتباط میان 
نقاط با و وبا در زمان له جاذب توسط تابع ۲ 
به‌صورت زير تقریب زده می‌شود (کوجاک و همکاران 
۰۷( 
)۸( (0 2 ربکا 

در روش تقریب موضعی پیش‌بینی» فرض می‌شود 
که تغییر ,26 بازمان در جاذب. در نقطه 4 با نقاط 
نزدیک آنها ( 2 ی,2ب1ح تم 25) یکسان باشد. 

در اینن‌جا م26 توسط ترتیب 0 ام از تابع 


11-1 1۳1-1 
‌ رک وت کی و » ‌ یط 2 و 5 وب 
20 1 


کا و1 
1,20 
11-1 
و ی ما ۳ ت 
ماع و 
ماو 
1,20 


)۹( 
بااستفاده از 9 تعداد ,جک و ِِِ برای مقادیر از 


پیش تعبین‌شده. ضرایب ] توسط معادلة زیر تعیین 


 < ۸۲ (۱۰(‏ 
که در این رابطه: 
)۱۱( ۱ و کدی )6 
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۹۸ 


۹9 


(رسم .نها و موول ری 1 س_ شوگ ول)<؟ 


و ۸ م‌اتریس ژاکوبین ۱/۲۳۱۵( +2<)۳ 


می‌باشد که به‌صورت زیر تعریف می‌شود: 


0 2 
و 2 با مه 4 4 
0 1 بو م1 1 
8 2 
بک- ر نع 2۶ ۶ ... 2 4 
۳۹ ۳ ۳۹ م1 ۳ - ۸ 
0 2 
بک مه 4 2 4 4 
رس و 1 رو ول تقو 1 


(۳ 


برای دست‌یابی به یک راه حل بهینه و مناسب. 
تاد سظر ماهس زب تایه بو سته ا ور 
زیر تبعیث کند [26] 


(0۱ (( ۳ < و 
0۱( 


شبکه‌های عصبی مصنوعی. شبکه‌های عصبی مصنوعی 
نگاشت يا تبدیل از یک فضای چندبعدی به فضای 
چندبعدی دیگر را انجام می‌دهند و قابلیت‌هایی از قبیل 
تاش اکن لیم ای عم اف 
توا زان دشن و کال دازا شی پات 
[۰3 ۲ ]. 

شبکه‌های چندلاية پیش‌خور یکی از مهم‌ترین 
ساختارهای شبکه‌های عصبی مصنوعی می‌باشند. به‌طور 
مقمون یماسا ادا شم هه ا مان راشتههاء: 
یس تفای تایه سوه کب تفا شا لاد 
ورودی. یک پا چندلابة پنهان و یک لاية خروجی 
می‌باشند. سیگنال ورودی در خلال شبکه و در مسیری 
روبه‌جلو به‌صورت لایه‌به‌لایه منتشر می‌شود. این نوع 
شبکه معمولاً با عنوان پرسپترون چند لایه(001۳) 
نامیده می‌شود. از میان روش‌های مختلف آموزش به 
روش پس‌انتشار خطاء الگوریتم لونبرگ- مارکوآد 
به‌دلیل هم‌گرایی سریع‌تر در آموزش شبکه‌های با اندازه 
فیح شاماد بر یی ای تسا شاه 


افتتتا: 
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ارزیاب ی کارایی روش پیش‌بینی موضعی ناشی از بازسازی فضای فاز و ... 


معیارهای ارزیابی مدل. معیارهای عددی که در تعیین 


می‌رود شامل: ضریب همبستگی (0 و جذر میانگین 


(۱۵ 2 ۳۳ 5 1 سم 
۹ ما 155 
۱ کیب 


که در این روابط مکی دادهٌ مشاهداتی و زمکی داده 
ماک نز مرف رایع 
هم‌چنین ی به‌ترتیب متوسط دادهٌ مشاهداتی 
و محاسباتی است و لا برابر با تعداد داده‌همای 
صحت سنجی می‌باشد. به‌طور کلی مقادیر ضصریب 
همبستگی نزدیک به عدد واحد و مجذور میانگین 
مربعات خطای نزدیک به صفر مبین عملکرد مناسب 


نتایج و بحث 

به‌منظور بررسی وجود آشوب در سری زمانی رواناب 
ابتدا زمان تأخیر هر یک از مقیاس‌های روزانه. هفتگی و 
ماهانه بااستفاده از روش میانگین اطلاعات متقابل 
(۸0۷) به‌دست آمده که در این روش اولین کمينة 
محلی از میانگین اطلاعات متقابل به‌عنوان زمان تأخیر 
درنظر گرفته می‌شود و مقدار بعد محاط با استفاده از 
روش نزدیک‌ترین همسایگی کاذب بااستفاده از 
زیربرنامه‌های نرم‌افزار تیسن محاسبه گردیده است. 
نمودار تغییرات اطلاعات متقابل سری زمانی رواناب در 
برابر زمان‌های تأخیر متفاوت در مقیاس‌های زمانی 
مختلف رسم گردیده است (شکل ۳. 

هم‌چنین مقادیر درصد خطای نزدیک‌ترین 
همسایگی کاذب به‌ازای مقادیر مختلف بعد محاط برای 


سال بیست و هشتم شمارهٌ دو ۱۳۹۵ 


محمد ذونعمت کرمانی- تحاطره امیرنحانی ۹۹ 


هر یک از مقیاس‌های روزانه هفتگی و ماهانه رسم و رواب روزانه سب 3 
0 ی 23۳ ۳4 
به‌عنوان بعد محاط مناسب انتخاب گردیده است ( شکل تس 
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شکل ۶ مقادیر نزدیک‌ترین همسایگی کاذب به‌ازای ابعاد محاط 
مختلف در سری زمانی رواناب 
رواناب ماهانه 3 
56 ام 
کم 5 ترسیم داده‌های سری زمانی داخل فضای فاز امکان 
۰ 04 7 نمایش سری‌های زمانی را در یک حالت مکانی می‌دهد. 
[ ۱ ۱ 
وم ۳۶ 4 درنهایت با توجه به مقادیر زمان تاخیر به‌دست امده از 
032 ت 
۱ 0 ِ روش تابع اطلاعات متقابل در شکل (۵) نمودارهای 
> ی ۱ ۳ 
0 ته فضای فاز دوبعد زما انات در مقباس ها 
‌ 3 ئّ 2 یار دوبدای سری: زمای زوابیاب دز معساسن ها ی 
زمان تاخیر (ماه) زمانی مختلف ترسیم شده اننحتا: همان‌طور که در شکل 


مشاه ی قرب اف ایک مقیاس مان از روتازی نا 


زماد اثات 
۱ هم‌چنین (۳)6 رواناب در زمان ] و (۳)۲ رواناب در 


زمان 1+0 و 7 مقدار زمان تن می‌باشد. 


سال بیست و هشتم شماره دو. ۱۳۹۵ نشریة مهندسی عمران فردوسی 


رواناب روزانه 
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شکل ۵ فضای فاز دوبعدی سری زمانی رواناب در مقیاس‌های 


زمانی روزانه» هفتگی و ماهانه (دبی برحسب مترمکعب بر ثانیه) 


رولاب روزائه 
>( سس 1 
7-220 روج«( 0 
2-40 5 4 ح 
ِ 5 
0260 سس / 8 2 ۹ 
0 / ۳ 3 
7-100 4 
20( 
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نشریة مهندسی عمران فردوسی 


ارزیاب ی کارایی روش پیش‌بینی موضعی ناشی از بازسازی فضای فاز و ... 


۲,09 )0( 
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1020( 


رواناب ماهانه 
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شکل ٩‏ نمودار مربوط به ارتباط میان تابع همبستگی( )0 و شعاع 
7 با افزایش بعد محاط (0) در سری زمانی رواناب 


در این مرحله به‌منظور محاسبهٌ تابع همبستگی 
رواتاب در عشقیاس‌های رسای شلف اشتدا نم فاز 
تغییرات (10801 در مقابل(108)7 ترسیم قیی کنو فقو 
نمودار شکل ۷ بش مسطحی در محدودهٌ (10801 
مشاهده می‌شود که در آن نسبت (1080(/086 به 
فاد تانق ماهتا این تایه این ک رس 
می‌توان در این محدوده انتخاب و مقدار بعد همبستگی 
موجود را محاسبه نمود. در ادامه به‌منظور تعیین بعد 
همبستگی نمودار توان همبستگی رسم گردید (شکل 
۹۱ 

همان‌طور که در شکل (۷) مشاهده می‌شود مقادیر 
توان همبستگی به‌ازای سری زمانی رواناب در هر سه 
مقیاس با افزایش بعد محاط زیاد شده و در نقاطی 
خاص به مقدار اشباع رسیده است. اشسباع تسوان 
همبستگی در مقدار معینی از بعد محاط بیانگر امکان 
وجود دینامیک قطعی و رفتار آشسوبناک در سیستم 
می‌باشد. بنابراین می‌توان نتیجه گرفت که سری زمانی 
رواناب احتمالاً دارای رفتار آشوبناک می‌باشد. 


سال بیست و هشتم. شماره دو ۱۳۹۵ 
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شکل ۷ نمودار تغییرات توان همبستگی با افزایش بعد محاط سری 


زمانی رواناب در مقیاس‌های زمانی مختلف 


نتایج سینت له از شکل (۷) در تعیین بعد 
همبستگی ۲ برای رواناب روزانه و ۶ برای رواناب 
هفتگی و ماهانهه همسو با مطالعات پیشین می‌باشد 
[9 ,]در اداته پررشتین الاک سر ژهانین 
بااستفاده از بزرگ‌ترین نمای لیاپانف صورت پذیرفت. 
شکل (۸) منحنی‌های پارامتر فاکتور کشش (5) درمقابل 
تعداد نقاط () را نشان می‌دهد» مقدار بزرگ‌ترین نمای 
لیاپانوف (۸) از برازش خحط بر روی این منحنی‌ها 
براساس روش حداقل مربعات و تعیین شیب آن صورت 


می‌پذیرد. 


سال بیست و هشتم. شماره دو ۱۳۹۵ 


رواناب روزانه 
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شکل ۸ نمودار مربوط به ماکزیمم نمای لیاپانوف سری زمانی 


رواناب در مقیاس‌های زمانی روزانه. هفتگی و ماهانه 


در صورت مثبت بودن بزرگ‌ترین نمای لیاپانوف 
در سری زمانی مورد بررسی. آشوبناکی آن سری زمانی 
محتمل می‌باشد. افزون بر این از ماکزیمم نمای لیاپانوف 
بر عم بر اوه اف نی سر ای زسانی اسفاده 
می‌شود. مقادیر افق پیش‌بینی( ۸) برای سری‌های زمانی 
رواناب در مقیاس‌های زمانی مختلف باتوجه به رابطءهٌ 
+2 ۸ محاسبه شده و نتایج آن در جداول (۳) ارائه 


گردیده انتختتا: 


نشریةٌ مهندسی عمران فردوسی 


۱۲ 


حدول ۳ مقادیر افق پیش‌بینی در مقیاس‌های زمانی مختلف سری 


زمانی رواناب 
رو 
بزرگ‌ترین نمای لیاپانف | افق پیش‌بینی 
روزانه 2 ۱9۸ 
هفتگی ۰۳۵ ۸۰ 
ماهانه ۰-۱۳5۲ ۳۹ 


همان‌گونه که از نتایج جدول (۳) برمی‌آید. مقدار 
مثبت بزرگ‌ترین نمای لیاپانوف در هر سه مقیاس 
تأییدی بر احتمال وجود رفتار آشوبناک در سری زمانی 
رواناب می‌باشد. 

پس از بررسی‌های صورت گرفته و مشاهده رفتار 
آشوبناک در سری زمانی رواناب حال می‌توان به بررسی 
روش‌های پیش‌بینی موجود به‌منظور دست‌یابی به 
روشی مناسب برای پیش‌بینی این سری زمانی پرداخشت 
که در تحقیق حاضر روش پیش‌بینی موضعی بدین 
منظور مورد استفاده قرار گرفته ‏ است. 

کلیةٌ مراحل محاسبات پیش‌بینی برای سال آخر 
از سری داده‌های ۱۲ساله. به‌عنوان دوره آزمون انجام 
گرفته است. به‌منظور فرآیند پیش‌بینی از مقادیر بعد 
محاط متخلف اسففاده گردیده است» تا بهترین پیش‌بیش 
در بعد محاط مناسب به‌دست آید. در جدول (4) مقادیر 
ضریب همبستگی (۲) و جذر میانگین مربعات خطای 
پیش‌بینی (1815۳) و بعد محاطی که بهترین پیش‌بینی 
به‌ازای آن صورت گرفته برای سری زمانی رواناب در 
مقیاس‌های زمانی روزانه» هفتگی و ماهانه ارائه شده 
است. هم‌چنین نمودارهای مقایسٌ مقادیر مشاهداتی در 
مقابل مقادیر پیش‌بینی شده نیز در شکل )٩(‏ نشان 


داده‌شده است. 


نشریة مهندسی عمران فردوسی 


ارزیاب ی کارایی روش پیش‌بینی موضعی ناشی از بازسازی فضای فاز و ... 


جدول ۶ ملفه‌های آماری مربوط به پیش‌بینی موضعی سری زمانی 
رواناب با ابعاد محاط مختلف در مقیاس‌های زمانی مختلف 
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شکل ٩‏ نمودار مقادیر مشاهداتی در مقابل بهترین مقادیر پیش‌بینی 
شده با روش موضعی در سری زمانی رواناب در ایستگاه پل کهنه, 
مربوط به داد سال آبی ۱۳۷۸-۱۳۷۷ 


سال بیست و ّ هستم شماره دوه ۱۳۹۵ 


محمد ذونعمت کرمانی- حاطره امیرحانی 


به‌منظور پیش‌بینی با شبکةٌ عصبی مصنوعی» شبکه 
یی و عون لا بابک لاید یی لکتورنیش 
آموزش لونبرگ-مارکوآد مورد استفاده قرار گرفت. 
به‌منظور مدل‌سازی. سال آخر به‌عنوان دادة آزمون و 
مابقی سال‌ها به‌عنوان داد آموزش سیستم انتخاب 
گردید. به‌منظور تعیین بهترین معماری شبکه. محاسبات 
در چندین حالت متفاوت باتوجه به روند افزایشی تعداد 
نرون‌های لاه میانی صورت پذیرفت و بابررسی 
میانگین مربعات خطاء بهترین ساختار شبکه تعیین 


شد. در ادامه ضریب همبستگی (۲) و جذر میانگین 
میات شا( 9۲ ور صتوی زره وش ذارستاض 
ای فاد و الا هدانی کو مقانا فاد سا ش شتتم 
نیز در شکل (۱۰) ارائه گردیده تاه 

درانتها نتایج حاصل از پیش‌بینی بادو روش 
موضعی و شبک؛ٌ عصبی در جدول (1) مورد مقایسه قرار 


گرفته است. 


جدول ۵ پارامترهای آماری دقت مدل‌های حاصل از شبکة عصبی مصنوعی در پیش‌بینی سری زمانی رواناب در مقیاس 


روزانه. هفتگی و ماهانه 
وس 2062 « | تعداد نرون لایذ مخفی الگوهای ورودی مقیاس زمانی 
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۲۳۱/۹3 ۰:۳ 31 (3- 000 ,(2- )00 ,(1- )0 


جدول 1 پارامترهای آماری مربوط به بهترین پیش‌بینی بااستفاده از روش موضعی و شبکهة عصبی مصنوعی. سری زمانی رواناب در 
تاش ها وتان مان 


سال بیست و هشتم شماره دو ۱۳۹۵ 


۳ روش موضعی روش شبکه عصبی 
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نشريهُ مهندسی عمران فردوسی 
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رواناب ملهانه 


ارزیاب ی کارایی روش پیش‌بینی موضعی ناشی از بازسازی فضای فاز و ... 


روائاب(میلی متر) 
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شکل ۱۰ نمودار مقادیر مشاهداتی در مقابل بهترین مقادیر پیش‌بینی به‌دست آمده از روش شبکة عصبی در سری زمانی رواناب. در ایستگاه 
پل کهنه. مربوط به دادهُ سال آبی ۱۳۷۸-۱۳۷۷ 


باتوجه به نتایج ارائه‌شده در جدول (1) مشخص 
می‌شود که روش پیش‌بینی موضعی توانایی بیشتری در 
پیش‌بینی سری زمانی رواناب در مقیاس‌های روزانه و 
هفتگی (۲-۰/۹۶ و ۰/۷۳-<1۷]9۳ برای روزانه و 
۵ و 3۷]5۲-۸۷/۷ برای هفتگی) نسبت به شبکهٌ 
عصبی مصنوعی (۲-۰/۹۲ و 2۰/۷۹ 15۳ برای 
مقیاس روزانه و ۲2۰/۹۶ و 2۸/۰۳ 80۷5 برای مقیاس 
هفتگی) داشته است. این در حالی است که شبکة عصبی 
مصنوعی (۲2۰/۷ و 2۳۰/۶۸ :10۷]5) عملکرد بهتری 
را نسبت به پیش‌بینی در مقیاس ماهانه نسبت به روش 
پیش‌بینی موضعی (۲-۰/۳۸ و 2۷۱/۳۸ :151) نشان 
می‌دهد. علت این امر را می‌توان در عدم آشوبناک بودن 
سری زمانی رواناب در مقیاس ماهانه دانست چرا که 
عملکرد این روش در مقیاس بزرگ‌تر نسبت به مقیاس- 
های زمانی کوچک‌تر به‌شدت افت پیدا کرده است. 
اهتمام به بحث عدم آشوبناکی مقیاس ماهانه بااستفاده از 


روش پیش‌بینی موضعی از این بابت ضروری است که 


نشریة مهندسی عمران فردوسی 


پیشتر وجود بعد همبستگی و توان مثبت لیاپانوف نوید 
احتمال آشوبناک بودن سری زمانی را داده بودند. علت 
این تناقض را می‌توان در عدم قطعیت روش‌های تعیین 
خحصوصیات دینامیکی یک سیستم قلمداد نمود که به این 
مسئله در بسیاری از مراجع اشاره شده است [5,11]. به 
همین دلیل حصول اطمینان از وضعیت آشوبناکی یک 
سری زمانی نیازمند بررسی و کنکاش در شیوه‌های 
متفاوت بررسی دینامیک‌های یک سیستم است. 
به‌طوری‌که عدم احراز مثبت بودن آشوبی بودن یک 
سری زمانی در هر یک از روش‌های متداول» آشوبناک 
بودن آن سری را زیر سژال می‌برد. افزون بر دلایل فوق 
یکی دیگر از علل بالقوة عدم احراز شرایط آشسوبناکی 
رواناب در مقیاس ماهانه را بااستفاده از روش پیش‌بینی 
موضعی می‌توان کمتر بودن تعداد داده‌های این سری 
نسبت به مقیاس روزانه و هفتگی برشمرد (جدول ۲ 
چرا که یکی از مهم‌ترین شرایط بررسی وضعیت 
آشوبناکی یک سیستم دینامیک» تعداد داده‌های زیاد آن 
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می‌باشد. علاوه بر اين لازم به ذکر است که مدل شبکة 
عصبی ساختاری متفاوت (جعبهٌ سیاه) را برای فرآیند 
آشوبناکی داده‌ها نیست و عملکرد مناسب‌تری را در 


تشر امش با شا ماهنت از رههان دنه است: 


نتیجه گیری 
در تحقیق حاضر به مدل‌سازی و پیش‌بینی سری زمانی 
رواناب در ایستگاه پل‌کهنه حوض آبریز قرهسو 
بااستفاده از روش‌های مبتنی بر تئوری آشوب و شبکهة 
عصبی مصنوعی پرسپترون پرداخته شد. در ابتدا 
آشوبناکی سری زمانی رواناب مورد بررسی قرار گرفشت 
و باتوجه به نتایم حاصل, امکان وجود آشوب بااستفاده 
از هر دو روش بعد همبستگی (مابین ۲ تا۶) و نمای 
لیاپانف (وجود نمای مثبت) در مقیاس‌های روزانه 
هفتگی و ماهانه مشاهده شد. هم‌چنین مدل‌سازی 
رواناب بااس‌تفاده از نظریٌ آشسوب و شبکه‌های 


منابع 


عصبی مصنوعی انجام شد. به‌منظور پیش‌بینی با 
به‌کارگیری مفاهیم تثوری آشوب. از روش پیش‌بینی 
موضعی که متکی بر بازسازی فضای فاز سیستم دینامیک 
است بهره‌گیری شد. به‌منظور بازسازی فضای فاز سری 
زمانی روش تاکنز مبتنی بر محاسبة زمان تأخیر مناسب و 
بعد محاط مورد استفاده قرار گرفت. نتایج حاصل از 
پیش‌بینی‌های صورت گرفته برتری روش پیش‌بینی 
موضعی را نسبت به روش شبکهٌ عصبی مصنوعی در 
مقیاس‌های روزانه و هفتگی بیان می‌دارد که این امر 
مژید وجود آشوب در مقیاس‌های فوق‌الذکر می‌باشد. 
شایان ذکر است که مدل شبکه؛ عصبی در پیش‌بینی 
رواناب ماهانه نسبت به روش پیش‌بینی موضعی توافق 
فاشفا ات کهدان اس ,هلت ماهتت ماوت اس و 
مدل غیرخطی در شبیه‌سازی سری زمانی می‌باشد. 
هم‌چنین اين مسئله نشانه‌هایی را از عدم آشوبناکی 
تسار سای اسان و سا مس ده 


. قبادیان ر. قربانی م.ع» خلج م. "بررسی عملکرد روش برنامه‌ریزی بیان ژن در روندیابی سبلاب رودخانة زنگمار درمقایسه با روش 
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